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ABSTRACT 

We study the distribution of gas pressure and entropy in eight groups of galaxies belonging to the 
ROSAT All-Sky Survey / Center for Astrophysics Loose Systems (RASSCALS). We use archival and 
proprietary XMM-Newton observations, supplementing the X-ray data with redshifts derived from the 
literature; we also list 127 new redshifts measured with the Gemini North telescope. The groups are 
morphologically heterogeneous in both the optical and the X-ray, and several suffer from superpositions 
with background galaxies or clusters of galaxies. Nevertheless, they show remarkable self-similarity in 
their azimuthally averaged entropy and temperature profiles. The entropy increases with radius; the 
behavior of the entropy profiles is consistent with an increasing broken power law with inner and outer 
slope 0.921q q5 and 0.421q q4 (68% confidence), respectively. There is no evidence of a central, isentropic 
core, and the entropy distribution in most of the groups is flatter at large radii than in the inner region, 
challenging earlier reports as well as theoretical models predicting large isentropic cores or asymptotic 
slopes of 1.1 as r ^ oo. The pressure profiles are consistent with a self-similar decreasing broken power 
law in radius; the inner and outer slopes are — O.TSIq'qs and — l.T^g g, respectively. The results suggest 
that the larger scatter in the entropy distribution reflects the varied gasdynamical histories of the groups; 
the regularity and self-similarity of the pressure profiles is a sign of a similarity in the underlying dark 
matter distributions. 



determination of elemental abundances. The cD galaxies 
in groups make up a larger fraction of the total baryonic 
mass of the system than do the cD galaxies in clusters 
(Helsdon et al. 2001; Lin & Mohr 2004); hence the prop- 
erties of the gaseous medium are more directly affected 
by the formation and evolution of the central galaxy. Fi- 
nally, because the typical velocities of galaxies in groups 
are comparable to the galaxies' internal (stellar) velocity 
dispersions, galaxy interactions through the merging insta- 
bility are common and very effective (Diaferio et al. 1993; 
Ramella et al. 1994). 

Recent studies suggest that in the X-ray, the distinc- 
tion between groups and clusters goes beyond an arbitrary 
temperature or mass boundary. Differences in the physical 
state of the intragroup medium make it difficult to view 
the < 2 keV systems simply as smaller, "rescaled" ver- 
sions of the more massive clusters. Of particular interest 
are the heating of the ICM by nongravitational processes 
such as supernova explosions, stellar winds, AGN activity, 
and shocks resulting from interaction with the surround- 

^ Based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA 
Member States and the USA (NASA). The XMM-Newton project is supported by the Bundesministerium fiir Bildung und Forschung/Deutsches 
Zentrum fiir Luft- und Raumfahrt (BMFT/DLR), the Max-Planck Society and the Heidenhain-Stiftung, and also by PPARC, CEA, ONES, 
and ASI. Also based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United 
States), the Particle Physics and Astronomy Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), 
the Australian Research Council (Australia), CNPq (Brazil), and CONICET (Argentina). 
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1. INTRODUCTION 

Groups are the intermediate structures between galax- 
ies and clusters of galaxies, and thus are crucial to an un- 
derstanding of how galaxies evolve in dense environments 
(Zwicky 1937; Zwicky & Humason 1960; Hickson 1982; 
Geller & Huchra 1983; RameUa, Geller, & Huchra 1989; 
White et al. 1999). In the X-ray, groups and clusters 
are customarily distinguished by the temperature of the 
hot, gaseous medium that makes up « 10% of their total 
mass (Ebeling, Voges, & Bohringer 1994; Ponman et al. 
1996; Mahdavi et al. 2000; Mulchaey 2000). Clusters, 
with temperatures > 2 keV, contribute only about 3% to 
the matter density of the Universe; but groups, with X-ray 
temperatures in the range 0.5-2.0 keV (corresponding to a 
mass range 10^^-10^'* Mq), contribute about twice as much 
(Ikebe et al. 2002; Reiprich & Bohringer 2002). Because 
of the cooler intracluster medium (ICM) temperature of 
the groups, the set of observable X-ray spectral lines is 
richer than it is in clusters, and allows a more accurate 
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ing large scale structure (Ponman, Cannon, & Navarro 
1999; Loewenstein 2001; Tornatore et al. 2003). Some- 
times called "preheating," these processes can leave a dis- 
tinct mark on the entropy distribution of the gas in groups, 
to a degree not observable in more massive clusters. For 
this reason groups provide a fossil record of this energy 
production during cosmic structure formation and galaxy 
evolution. 

In this study we focus on the distribution of entropy and 
pressure in cnglit systems drawn from the RASSCALS sur- 
vey of nearby galaxy groups. In §2 and §3 we discuss the 
optical and X-ray data we have gathered for this study, as 
well as the data reduction procedures used in our analysis. 
In §4 we describe our attempts to fit self-similar profiles 
to the entropy and the pressure distribution. In §5 we 
describe each system in detail. In §6 we summarize our 
conclusions. 

2. OPTICAL DATA 

2.1. Sample Selection 

At the focus of our study are groups with emission- 
weighted temperature ^ 0.5 — 2 keV. The X-ray 
luminosity-temperature relation (Ikebe et al. 2002) sug- 
gests that these groups should have Lx < 10^^ erg s~^ in 
the 0.1-2.4 keV energy range. ^ 

To build our sample, we begin with the ROSAT All-Sky 
Survey-Center for Astrophysics Loose Systems (RASS- 
CALS), a statistically complete, magnitude-limited cat- 
alog of optically identified groups with 0.01 < z < 0.04 
(Mahdavi et al. 2000). This catalog contains 260 groups, 
of which 43 have X-ray emission detected in the RASS. 
To conduct a more detailed study of the properties of 
these groups, we select X-ray emitting RASSCALS with 
Lx < 10*^ erg s~^ that have been observed with the 
XMM-Newton observatory. The properties of the sample 
are shown in Table 1. 

The membership of the original RASSCALS groups was 
established spectroscopically with a completeness limit 
ruR « 14.4 (Ramella et al. 2002). All but one of the 
groups we have selected for analysis have also been the 
targets of deeper redsliift surveys, by Mahdavi & Geller 
(2004) and Rines et al. (2003) (complete to ttir « 15.4), 
by Zabludoff & Mulchaey (1998) (70%~95% complete to 
different magnitude limits depending on the group) , or by 
Pinkney et al. (1993) (imknown completeness). None of 
the groups had significant overlap with the current release 
of the Sloan Digital Sky Survey (Abazajian et al. 2004). 
In addition, we have performed optical observations which 
we describe below. 



2.2. New Optical Spectroscopy 
We observed three groups, RGH 80, HCG 97, 



and 



NRGbl84 with the 8m Gemini North telescope on Mauna 

Kea, Hawaii. The Gemini multi-object spectrograph, 
CMOS, was used to measure rcdsliifts for galaxies as faint 
as mji w 20 within 6.5' of the center of each system. 

The first step was the selection of targets for spec- 
troscopy. We used GMOS in queue imaging mode to ob- 
tain a 2 X 2, 13' X 13' i?-band mosaic around the central 



galaxy in each group; the length of each imaging expo- 
sure was 10 minutes. The SExtractor package (Bertin & 
Arnouts 1996) then separated galaxies from stars, and cal- 
culated the magnitudes and half-light radii of the galax- 
ies. It would have taken an inordinate amount of time to 
measure redshifts for all mn < 20 galaxies in each image. 
Therefore, galaxies were sorted by their half-light radii, 
with the assumption that regardless of magnitude, galax- 
ies with larger half-light radii would be more likely to have 
low redshifts and therefore belong to the group. All three 
groups are at a small enough redshift that foreground con- 
tamination was not an issue with this selection procedure. 

In each piece of the mosaic, we selected the 20 galax- 
ies with the largest half-light radii and without previ- 
ously measured redshifts. Thus there were a total of 80 
targets for NRGbl84 and HCG 97. For RGH80, we in- 
cluded an additional central mask, for a total of 100 tar- 
get galaxies. We designed one slit mask per image to 
be used with the GMOS B600 grating; the slit size was 
0.5"(RGH80 and HCG 97) or 0.75"(NRGbl84) achieving 
a resolution of 4 — 6A at 4000A. The wavelength coverage 
was 4000 - 6OOOA, with a shift of - lOOOl in either direc- 
tion depending on the position of the galaxy on the focal 
plane. We used two 30 minute exposures per slit mask 
(a total of one hour), taking arc lamp exposures between 
each set of two exposures to obtain accurate wavelength 
calibrations. Because absolute photometric calibration is 
not required for redshift measurements, we did not apply 
flat-fielding or fiux calibrations to the spectra. 

The Gemini GMOS package for IRAF^ was used to cal- 
culate the wavelength solutions and to reduce the multi- 
object observations into one-dimensional spectra. The 
RVSAO package (Kurtz & Mink 1998), incorporating the 
methods of Tonry & Davis (1979), allowed us to measure 
redshifts by maximizing the cross-correlation of the spec- 
tra with absorption- and emission- line templates. The es- 
timation of errors in the optical velocities is described in 
detail by Kurtz & Mink (1998), who determine the rela- 
tionship between the shape of the cross-correlation func- 
tion peak and the 68% velocity confidence interval using 
galaxies with known velocities. Our results are in Table 2. 

2.3. Membership 

We assemble a galaxy catalog by combining all the avail- 
able redshift surveys in the direction of each group, and 
restricting ourselves to galaxies within 2 Mpc of the group 
center. The mean redshift z and velocity dispersion a of 
each group are defined as follows (Danese, de Zotti, & di 
TuUio 1980): 
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(A^-l)(l + z)2 

The subtraction of the root-mean-scjuarc velocity error in 
the second term is an attempt to remove the contribution 
of the measurement uncertainties ej to a. In practice, this 
correction may be inaccurate when the underlying veloc- 



ity distribution is nongaussian, the sample size is small 
^ We assume i?o = 70 km Mpc~^, Oq = 0.3, and Oa = 0.7 throughout the paper. 

^ For a description of the Gemini data reduction pa<;kage see http://www.gemini.edu/sciops/data/dataSoftwaje.htmI. 
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{N < 20), and the velocity dispersion is of the same or- 
der of magnitude as e,. The inaccuracy is not a source 
of concern for our sample, but may be relevant for some 
extremely low-velocity dispersion groups with cr w 70 km 
(Mahdavi ct al. 1999, 2000). 
To determine group membership, we then use the 
"sigma-clipping" (Zabludoff, Huchra, & Geller 1990; Mah- 
davi & Geller 2004) algorithm to reject outliers, or galaxies 
unlikely to be bound to the group. This algorithm consists 
entirely of making sure that no group member is separated 
from its nearest neighbor in velocity space by more than 
the velocity dispersion a of the group. 

3. XMM-NEWTON OBSERVATIONS 

We conducted XMM-Newton Observations of three of 
the eight groups, while the remaining five arc publicly 
available in the XMM-Newton Data Archive. Table 3 con- 
tains the details of the observations. 

The initial steps of the data reduction are similar to 
the procedure described in Zhang et al. (2004) and 
Finoguenov et al. (2003). The first important aspect is 
the removal of flares, which can significantly enhance the 
detector background, severely limiting the detection of low 
surface brightness features. Thus, for the group analysis, 
using flare free observing periods is critical. At energies 
above 10 keV the particle background dominates the de- 
tected counts; there is little X-ray emission from our tar- 
gets because (1) the telescope efflciency is quite low at 
energies > 10 keV, and (2) the temperature of the objects 
in our study is less than 2 keV. We use the 10 15 kcV 
energy band (binned in 100s intervals) to monitor the par- 
ticle background and to excise periods of high particle flux. 
In this screening process we use the settings FLAG=0 and 
PATTERN< 5 for the pn detector on XMM-Newton. We 
reject time intervals affected by flares by excising periods 
where the detector count rate exceeds the mean quiescent 
rate by more than 2(t. To produce the broad-band images 
and hardness ratio maps, we have also used the MOSl 
and M0S2 events with PATTERN< 13 and FLAG=0. To 
reduce the widths of the gaps in the pn broad band im- 
ages, we included photons near the pn-CCD borders, near 
bad pixels, and near offset columns. Data from the pn de- 
tector were used exclusively in the spectroscopic analysis 
described below. 

Because most of the observations analyzed here were 
performed using a short integration frame time for pn (Full 
Frame Mode), it is important to remove the out-of-time 
events (OOTE) for accurate imaging and spectral analy- 
sis. We used the standard product of the XMMSAS 5.4 
epchain task to produce the simulated OOTE flle for all 
the observations and scale it by the fraction of the OOTE 
expected for the frame exposure time, (Table 3). 

The cosmic component of the background consists of 
emission from the Galaxy as well as the extragalactic Cos- 
mic X-ray Background (CXB). Observations of blank fields 
also contain both components. Provided that the expo- 
sures are done with the same instrumental set-up (e.g. 
with a particular filter) the spectra of the CXB are the 
same for both the target and the blank field. The Galac- 
tic component acts as an absorber and emitter, thus it is 
important to choose a similar absorbing column density for 

* For further details of XMM-Newton processing we refer the reader 



both the target and background data sets, so that the ex- 
pected background spectra are similar. In addition, there 
are variations in the Galactic emission on small scales. 

The vignetting correction is performed taking into ac- 
count the source extent, recent vignetting calibration 
(Lumb et al. 2004), and the pn response matrices released 
under XMMSAS 6.0. The residual systematic uncertainty 
of the flux is below 4% for the pn (Lumb et al. 2004). 

The analysis of each group consists of two steps: (1) 
construction of separate wavelet-decomposed maps of es- 
timated surface brightness in the 0.5-2.0 keV band (/) and 
of emission- weighted temperature (T/), which may be used 
to derive the entropy and pressure "integrated" along the 
line of sight; and (2) use of the wavelet-decomposed maps 
to identify contiguous regions from which we extract X-ray 
spectra for independent model-fltting. 

The background-subtracted wavelet maps are based on 
photon images corrected for instrumental effects. The sur- 
face brightness map is constructed using the technique 
described in Vikhlinin et al. (1998). A hardness map 
of the emission then results from dividing the wavelet- 
reconstructed images in the 0.5-1 and 1-2 keV bands. 
The hardness ratio is a monotonic and unique function of 
the emission- weighted temperature, as long as the group 
redshift is well known (true for our sample). The advan- 
tages of using wavelets include the ability to remove ad- 
ditional background by spatial filtering and a control over 
the statistical significance of the detected structures. We 
use the "a trous" method of wavelet image reconstruction 
with scales from 8" to 4', applying a four sigma detection 
threshold and retaining emission to a l.Tcr detection limit. 
Occasionally, the wavelet algorithm generates small scale 
discontinuities in the reconstructed image. We remedy this 
effect by applying additional smoothing before producing 
the hardness ratio maps. 

We also construct maps of the entropy and pressure "in- 
tegrated" along the line of sight, using the surface bright- 
ness image and the emission-weighted temperature maps. 
The use of the word "integrated" is not strict, because 
these quantities do not have the dimensions of entropy or 
pressure times distance, and they serve as heuristic aids 
only. In studies of the intracluster medium, the common 
definitions of entropy and pressure are kT/n^^^ and nkT, 
respectively, where n is the gas density and kT is the tem- 
perature in units of keV (see §4 for further details). Be- 
cause the emission measure is proportional to n^, we use 
\/7 as an estimate of the gas density "integrated" along the 
line of sight. The "integrated" entropy is calculated simply 
as T//^^/"^, and the corresponding pressure as VJTj. The 
quantities shown in these maps are not simply related to 
the physical, unprojected entropy and pressure, and they 
suffer from a number of degeneracies, with metallicity- 
density being the strongest (a significant fraction of the 
group emission results from line emission). Nevertheless, 
the maps indicate the regions of primary interest for fur- 
ther spectroscopic analysis, in which most of the degen- 
eracies are removed. In general, we expect the entropy 
to increase monotonically with radius (Metzler & Evrard 
1994), while the pressure should decrease. We stress that 
the "integrated" pressure and entropy are used only as 
tools to guide further spectroscopic analysis. They are 

http: / /wave.xray.mpe.mpg.de/xmm/cookbook/general. 
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not used in the fits we discuss in §4. 

The second, spectroscopic part of the analysis uses a 
mask file, constructed from both hardness ratio and sur- 
face brightness analysis described above. The first appli- 
cation of this technique is in Finoguenov et al. (2004). In 
our analysis we select regions with similar spectral prop- 
erties. We combine the regions so that counting statistics 
are not the limiting factor in our derivation of the group 
properties. We use the wavelet-based maps to identify re- 
gions with similar X-ray colors and intensity levels. To 
generate the mask file for use in the further spectral anal- 
ysis, we sample the changes in the intensity and hardness 
ratio at the precision allowed by the statistics of the data. 
We then examine each of the isolated regions with approx- 
imately equal color and intensity, imposing the additional 
criterion that the reg ions should be larger than the PSF 
width (15") and contain more than 300 counts in the raw 
pn image. A sample mask file is shown in Figure 1. 

The spectral analysis was performed using single- 
temperature models and the APEC plasma code, fitting el- 
emental abundances from O to Ni as one group and assum- 
ing the photospheric solar abundance ratios of (Anders & 
Grevesse 1989). It should be noted that evidence for mul- 
tiphase gas in groups of galaxies exists (Buote et al. 2003), 
such that gas at each projected radius may vary in tem- 
perature by a factor of w 2; however, if the variations are 
due to a smooth radial trend in the plasma temperature 
(i.e., a temperature profile), a series of single-temperature 
fits across the group should roughly reproduce the general 
shape of the underlying pressure and entropy distributions. 

In our analysis, we paid special attention to the issue of 
background estimation. We employ a double background 
subtraction technique, following Zhang et al. (2004). We 
used the region with radius 12-16' and 2-12 keV band to 
estimate the quality of the background subtraction using 
blank fields. This region is free of significant X-ray emis- 
sion from the groups. We fix the shape of the residual 
background component, and also add a 0.2 keV thermal 
component (APEC with solar element abundance) to ac- 
count for a possible variation in the Galactic foreground, 
allowing the normalizations of both components to be fit. 

To estimate the true (as opposed to "integrated" ) pres- 
sure nkT and entropy kTn~^^^ in each region, we need 
measurements of both the gas density n and the temper- 
ature T. The spectral fit products, however, are kT and 
n^V, where V is the volume of the emitting region. We 
do not have an exact knowledge of V. Because the mask 
generation technique identifies interesting regions by their 
wavelet-decomposed properties — producing region maps 
like Figure 1 — a straightforward deprojection of the spec- 
tra is not possible. We therefore need to estimate the 
length of the column for each selected two-dimensional re- 
gion on the sky. We assume that the gas emits uniformly 
along the line of sight as shown in Figure 2. With these 
approximations, the longest length through each volume is 
L = 2-\/i?2 — Ri, and the mean three dimensional distance 
of the region from the cluster center is r = (i?i + i?2)/2. 
The volume of the region is then 2SL/3, where S is the 
area of the region. A similar procedure is used in Henry, 
Finoguenov, & Briel (2004). 




Fig. 1. — A sample extraction template for HCG 97. The spectra 
in each region are analyzed independently. White circles represent 
excised unrelated sources. In this particular image, darker colors 
represent higher emission-weighted plasma temperatures. 




Ri 

^ R2 

V 

Observer 

Fig. 2. — Simplification of the group geometry. The emitting re- 
gion is to the left of the dashed line, and has a minimum and max- 
imum projected distance Ri and R2, respectively, from the group 

center. The distance along the line of sight is then 2^ i?| — R^. 

Figures 3 and 4 show the gas density, temperature, en- 
tropy, and pressure distribution generated by the wavelet 
decomposition technique. The locations of the member 
galaxies also appear in these figures. Table 4 lists the ba- 
sic physical properties of the X-ray emitting medium. 

4. ENTROPY AND PRESSURE PROFILES 

The distribution of entropy in systems of galaxies has 
been used to argue for the injection of "extra energy" 
into the intragroup medium during the formation process. 
Several authors (Ponman et al. 1999; Lloyd-Davies, Pon- 
man, & Cannon 2000) have argued that in addition to 
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Flo. 3. — X-ray and optical images of the eight groups in our sample. (Top) Wavelet decomposition of the X-ray surface brightness in 
the 0.5-2.0 keV energy band. Each image is 19' X 19'; the grid cells are 200 kpe on a side. The small boxes show the area enlarged below. 
(Bottom) Palomar Observatory Sky Survey images of the central region of each group, with X-ray surface brightness contours superimposed 
from the above images. Each image is 5.2' x 5.2' in size; the grid cells are 50 kpc on a side. 
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Fig. 4. — {Top) Emission- weighted temperature distribution of the sample from wavelet analysis. The range in temperatures shown is 0.5 
keV (white) to 2.3 keV (black), with a linear scale. Member galaxies are shown as crosses. Each grid cell is 200 kpc on a side. (Bottom) 
Distribution of the entropy "integrated" along the sight, estimated from wavelet analysis. The range in integrated entropy shown is 10 (white) 
to 110 (black) in arbitrary units, with a linear scale. The contours show lines of constant gas pressure "integrated" along the line of sight; 
the innermost contours always show higher pressure, with each additional contour showing a decrement in pressure from the previous one of 
a factor of 1.5. 



7 



gravitational collapse, processes such as supernova explo- 
sions make a contribution to the thermodynamic state of 
the gas. This contribution takes the form of a "floor" in 
the inner entropy distribution of the intragroup medium 
as reported in these studies. In other words, gas with 
higher entropy than expected from pure gravitational col- 
lapse appears to be pooled at the center of groups ob- 
served by the ROSAT satellite. Ponman, Sanderson, & 
Finoguenov (2003) use a larger sample of 66 systems ob- 
served by ASCA to study this effect in greater detail. They 
find evidence of an entropy floor in these observations, 
with the notable exception that in lower temperature (less 
massive) systems the effect is much less pronounced than 
in higher temperature (more massive) ones. The entropy 
profiles of the less massive systems are better described by 
a simple power law. 

Further complicating the picture, Finoguenov et al. 
(2002), using ASCA, foimd that gas at rsoo (the radius 
within which the mean matter density is 501 times the crit- 
ical density of the universe) also exhibits excess entropy 
relative to the pure gravitational value. This outer en- 
tropy excess seems consistent with models in which shock 
heating and galactic winds are also a major contributor to 
the groups' dynamical evolution (Dos Santos & Dore 2002; 
Voit & Ponman 2003). In addition, Ponman et al. (2003) 
found that the entropy profiles, once scaled to T°-^^, are 
in agreement with each other, hinting at the possible uni- 
versality of the entropy enrichment phenomenon. Sonic; 
XMM-Newton observations are consistent with this result 
(Pratt & Arnaud 2002). Next we revisit all of these results. 

4.1. Calculation of the Profiles 

Figures 5 and 6 show the pressure and entropy as a func- 
tion of true (as opposed to projected) distance from the 
center for our sample. The true distances were estimated 
not via deprojection but using the method described in 
§3 above. We define the gas entropy similarly to the ear- 
lier discussions of the intragroup medium (Ponman et al. 
2003): 

5 = Tn-t, (3) 

where n is the particle number density and T is the tem- 
perature in keV. As a result, S has dimensions of keV 
cm^. The traditional, thermodynamical entropy is then 
logs' -H C, where C is an arbitrary constant. 

We use the spectral fits discussed in §3 to plot the en- 
tropy as a function of distance from the cluster center. The 
fits arc carried out in each region with the abundance, tem- 
perature, and gas density as free parameters, and with the 
hydrogen column fixed to the Galactic value from 21cm 
surveys (Dickey & Lockman 1990). The results are then 
used to calculate the entropy via equation (3). We add 
an additional 10% uncertainty to the entropy and pres- 
sure calculated in this manner to account for inaccuracies 
involved in our method of estimating distances from the 
center (§3); the volume estimates are more uncertain than 
the distances. 

Because the entropy and pressure profiles in the groups 
appear similar, we test whether simple empirical models 

can describe their shapes, and whether the same entropy 
and pressure profiles can describe all the groups in our 

^ Sherpa was used for fitting the radial profiles only; details on the s 



sample. Rather than scaling the data by other derived 
values such as the temperature and rsoo , we perform mul- 
tidimensional fits (in up to 18 dimensions) as required to 
compare the data with the self-similar models. 

4.2. Single Power Laws 

To begin, we test whether the entropy and pressure pro- 
files can be described by a single power law. From a the- 
oretical point of view, such a profile seems unlikely, be- 
cause N-body simulations (Navarro, Frenk, & White 1997; 
Ghigna et al. 2000) suggest that the mass distribution 
in dark haloes is not well fit by a single power law. The 
pressure and the entropy distrbution depend on the tem- 
perature, and the temperature profile in turn depends on 
the overall matter density distribution. Thus it would be 
surprising if the resulting entropy and pressure were scale- 
free. However, noise may well make the observations con- 
sistent with a single power law, and it is important to test 
for this possibility before considering more complicated 
models. 

We wish to discover whether the pressure and entropy 
distributions in our sample of groups are self-similar. 
Therefore, in fitting single power laws to the data, we con- 
strain the slope of the power laws to be the same for all 
eight groups: S{r) cx r^' for the entropy, and P{r) cx r^^ 
for the pressure. Each fit has 9 free parameters: one slope 
and eight normalizations. We minimize the statistic via 
the built-in implementation of simulated annealing (Press 
et al. 1992) in the Sherpa data analysis package^. Sim- 
ulated annealing is a compromise between downhill mini- 
mization (which is fast but may not always find a global 
minimum) and monte-carlo minimization (which is slow 
but has a greater likelihood of finding a global minimum). 

The results appear in the top halves of Figures 5 and 6. 
The fits are not acceptable, yielding reduced of 2.42 and 
3.32 for the entropy and pressure data, respectively. We 
conclude that neither the entropy nor the pressure profiles 
are consistent with a single universal power law. 

We can proceed either by treating the groups as indi- 
vidual entities — fitting independent single power laws to 
each system — or by considering more complicated univer- 
sal models. The first approach yields an acceptable fit for 
the entropy profiles of Abell 2634 (^^ = 0.21 ±0.10), NGC 
5129 ifis = 0.79 ±0.08), HCG 97 (/x^ = 0.59 ±0.04), RGH 
80 = 0.62 ±0.03), SS2bl53, (n, = 0.72 ±0.04), and 
SRGbll9 {fj,^ = 1.1±0.1); Abell 194 and NRGbl84 do not 
yield acceptable fits. For the pressure profiles, single power 
laws yield a good fit to Abell 2634 {fip = -1.2±0.1), NGC 
5129 {fip = -1.0 ± 0.1), Abefl 194 {fip = -0.92 ± 0.11), 
NRGbl84 = -0.73±0.05), SS2bl53 (np = 1.26±0.04), 
SRGblig inp = 1.1 ± 0.09); HCG 97 and RGH 80 are not 
consistent with a single power law. Interpreting these in- 
dividual slopes is by no means straightforward, especially 
in light of the fact that we may not be sufficiently sam- 
pling the central regions of the more distant systems. It is 
far more instructive to consider whether a more detailed 
universal profile can describe all the data. 

4.3. Broken Power Laws 
.re may be found at http://asc.haj-vaj:d.edu/sherpa/ 
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Fig. 5. — Pressure distributions for the eight groups in our sample, (top) Single-power law fit to the entire sample, with the slope contrained 
to be same for each group; the reduced is 3.32. (bottom) Broken power law fit; the reduced is 1.18. 
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Fig. 6. — Entropy distributions for the eight groups in our sample, {top) Single-power law fit to the entire sample, with the slope contrained 
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We test the self-similarity of the entropy and pressure 
profiles by fitting broken power laws of the form 



P{r) = Po 



(r/lGO kpc)"" r < Tp 

(rp/100kpc)"''-'>'f (r/100 kpc)''f r > rp 



Similarly, we fit entropy profiles of the form 



Sir) = So 



(r/lOOkpc)"' 
(rJlOOkpc)"-' 



'(r/lOOkpc)''- 



r < Ts 
r > r. 



(4) 



(5) 



For these self-similar models we conduct two separate 
analyses: a pressure fit and an entropy fit. The pressure 
fit has 18 free parameters: the slopes ap and jp common 
to all the groups, and eight different pairs {rp, Pq) for each 
group. Similarly, the entropy fit has its own set of 18 pa- 
rameters. 

Minimizing using the same technique as in the pre- 
vious section yields good fits; these are shown in the bot- 
tom half of Figures 5 and 6. In Table 5, we estimate the 
goodness-of-fit (the probability of observing > Xmin 
the model is correct) as the definite integral of the x dis- 
tribution, q = r(z//2,x^/2)/r(z//2), where r{a,x) is the 
incomplete gamma function, T{a) = r(a, 0) is the gamma 
function, and v is the number of degrees of freedom (i.e., 
the number of data points minus the number of free pa- 
rameters minus one). 

After finding the minimum, we calculate errors on the 
best-fit parameters by using likelihood ratio tests (Lupton 
1993). This method fully takes the correlations among 
the parameters into account. The test involves repeatedly 
comparing the overall minimum with x^ minimized as- 
suming many different fixed values for the parameter in 
question. The set of differences between the global min- 
imum and the constrained minima has a distribution 
with one degree of freedom asymptotically, and this fact 
may be used to estimate confidence intervals. For each pa- 
rameter we re-minimize 100 times around the best-fit 
value; the error estimation procedure thus involved 1800 
re-minimizations in all. 

4.4. Results 

We find that a broken power law provides a good de- 
scription of the combined pressure profiles of the eight 

groups (see Table 5 and Figure 5). The pressure bro- 
ken power law has a shallow slope — 0.78^q q3 near the 
group center, steepening to — I.T^q J at the outer edges; 
the goodncss-of-fit is an acceptable q = 0.15. The agree- 
ment is not perfect: NRGbl84 is entirely consistent with 
a single, rather than a broken power law, and SS2bl53 ex- 
hibits a somewhat steeper inner pressure profile than the 
other seven groups. Otherwise, the notion that the groups 
are self-similar systems differing only in a radius and ab- 
solute pressure scaling is supported by the XMM-Newton 
data. 

The entropy profiles (Figure 6), however, make the pic- 
ture more complicated. The broken power law fit to the 
combined data has fit quality q = 0.0095 ; the model 
is acceptable at 3a. Excluding the binary cluster Abell 
194 would bring the goodness-of-fit to a wholly acceptable 
q = 0.136, and would not affect the best-fit slopes sub- 
stantially. The inner power law slope, 0.92tQ;g5, is steeper 
than the outer entropy slope 0.42tQ;Q4, meaning there is no 



evidence for an constant-entropy "floor" as close as 10 kpc 
from the center of each group. Rather, in contrast to pre- 
vious observations with ROSAT and ASCA, the entropy 
profile shows a tendency to steepen towards the center of 
each group. There is an entropy deficit at larger radii rela- 
tive to the extrapolation one would obtain from the inner 
power law, a new result not seen in data from the older 
X-ray observatories. 

It is instructive to compare our data with theoretical 
models. Tozzi & Norman (2001) make detailed predictions 
for the entropy profiles of groups and clusters of galaxies 
in two different scenarios: with preheating (where a signif- 
icant initial entropy excess exists in the systems at t « 1 
Gyr), and without preheating (where the gas dynamics is 
chiefly determined by shocks and gravitational heating). 
In the first scenario, they find that more massive clus- 
ters develop isentropic cores within ~ 0.05r20O5 and that 
less massive groups develop relatively larger cores (within 
~ 0.5r2oo). By contrast, in the scenario without preheat- 
ing, groups and clusters both tend to develop power law 
entropy distributions, with S{r) oc r^'^. 

To test these scenarios, we combine our entropy profiles 
in several ways (Figure 7). In (a), we show all the groups 
for which the preferred fit is a broken power law — that is, 
the groups for which is neither a lower nor an upper 
limit. Here the radii are scaled by r^, and the entropies 
by S{rs)- As might be expected, this combination of pro- 
files exhibits little sc;atter around the best-fit line, because 
the highest contributors to x^— Abell 194 and SRGbll9 — 
are not included. In (b) we show the profiles for all eight 
groups, scaled by rsoo and (S), the mean entropy within 
^500 as listed in Table 4. Despite the difference in scal- 
ing, the data are still well described by our best-fit broken 
power law with fixed slopes = 0.92 and 7^ = 0.42 — in 
fact, the fit quality is higher, because the additional uncer- 
tainty introduced through division by (S) decreases the re- 
duced x^- We compare these profiles with the predictions 
of the Tozzi & Norman (2001) preheating model, where we 
assume that (5) = 5'(x,/2, where is the entropy of the 
last accreted gas shell. The model is rejected, because it 
predicts a large isentropic core that we do not observe and 
an outer profile that rises more quickly than ours. Vary- 
ing {S)/Sx, arbitrarily in either direction does not change 
this result. Finally, in (c) we show the eight groups again, 
with the radii scaled by rsoo, but with the entropies scaled 
by (r/keV)" *'^ a relation Ponman ct al. (2003) derive 
for their sample of 66 virialized groups and clusters. This 
scaling agrees with our Ug = 0.92 and 7s = 0.42 fits as 
well, but is clearly inconsistent with both the preheated 
and non-preheated Tozzi & Norman (2001) models — both 
models are rejected with q < 10~^. 

Thus there are two key differences between our results 
and the Tozzi & Norman (2001) models: (1) we do not 
observe any isentropic cores, and (2) we find an entropy 
decrement at radii > O.Srsoo with respect to the extrap- 
olation from data at radii < O.Srsoo. Neither of these 
results is unique to our sample: both effects are also clear 
in the high quality Chandra observation of the NGC 1550 
group (Sun et al. 2003); the lack of an isentropic core 
is also clear in the Chandra and XMM-Newton observa- 
tions of ESQ 306170 (Sun et al. 2004) and in a few of 
the groups studied by Ponman et al. (2003). These stud- 
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Fig. 7. — The entropy profiles combined using (o) our fit results, but leaving out groups where S{r) is consistent with a single power law; 
(b) all groups, with the radii scaled by rsoo and the entropies scaled by the mean entropy {S) within rsoo; and (c) all groups with the radii 
scaled by rsoo and the entropies combined according to the scaling relation derived by Ponman et al. (2003). The solid curve in all three 
figures shows the broken power law derived from our analysis, as = 0.92 and 73 = 0.42, with the breaking radius and normalization adjusted 
to fit the data in (6) and (c). The short-dashed line in (a) shows the continuation of the as = 0.92 power law. The long-dashed curve in (6) 
shows the prediction of the Tozzi & Norman (2001) preheating model. The dotted curve in (c) shows the (rejected) best-fit power law with 
slope fixed at 1.1. 
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ies employ deprojection techniques very different from our 
mask-based approacli; thus it is unhkely that we are wit- 
nessing a pecuUar effect due to the instrument or due to 
our data reduction technique. 

In the inner regions, the disagreement with older studies 
may be due to the improved spatial resolution of Chan- 
dra and XMM-Newton. In studies such as Ponman et al. 
(2003), the distribution of the gas in groups is assumed 
to follow a single /3-model distribution, with gas density 
Pg oc (l+r2/r2)-''''3/2 (Sanderson ct al. 2003). Higher res- 
olution data indicate, however, that the gas density rises 
more steeply towards the center due to the influence of the 
dominant central galaxy (Helsdon & Ponman 2000; Lewis, 
Buotc, & Stocke 2003; Buote et al. 2003; Sun et al. 2003, 
2004), with n better described by a broken power law or 
by two superposed /3-models. In this scenario, as r — > 0, 
the entropy S oc T/ri^/^ declines more rapidly than it does 
for an equivalent single /3 model, unless Vc is very small. 

The disagreement with the theoretical models predicting 
large isentropic cores is likely related to the assumed equa- 
tion of state of the gas. In such models (Tozzi, Scharf, & 
Norman 2000; Tozzi & Norman 2001), the initial gas den- 
sity profile follows a polytropic distribution, n oc p^/'' oc 
7"i/('?-i)^ with T] being the adiabatic index. Given this 
functional relationship, the temperature gradient dT/dr 
is inevitably negative everywhere, because the gas density 
must decline with radius. This scenario, when combined 
with an initial entropy excess (due; to supernova-driven 
galactic winds, for example), yields an isentropic core that 
grows as the group evolves. However, the temperature dis- 
tribution in this preheating scenario is in strong contrast 
to the data, where, at least within 0.3r5oo, the temper- 
ature typically exhibits a positive gradient, rising mono- 
tonically from the center of groups as well as relaxed rich 
clusters of galaxies, and declining or remaining constant 
only outside 0.3r5oo (Pratt & Arnaud 2002; Zhang ct al. 
2004). Thus there is a fundamental problem with the 
inner temperature structure predicted by the polytropic 
preheating models. In the era of ROSAT and ASCA, 
this disagreement was sometimes explained with the argu- 
ment that the inner regions of clusters and groups contain 
cooling flows, which disrupt the equilibrium physics, and 
which therefore need to be excised (Markevitch 1998) or 
modeled as two-temperature systems (Ettori 2000). How- 
ever, Chandra and XMM-Newton show that hydrostatic 
models — in many cases including the region within the 
cooling radius — can fit clusters with rising temperature 
profiles (Schmidt, Allen, & Fabian 2001; Allen, Schmidt, 
& Fabian 2002; Ettori et al. 2002; Buote et al. 2003; 
Zhang et al. 2004). 

Therefore, rising temperature profiles must be a part 
of any successful preheating model. The preheating mod- 
els of Muanwong et al. (2002), Tornatore et al. (2003), 
and McCarthy et al. (2004), for example, do away with 
the polytropic assumption, and therefore can reproduce 
the positive inner temperature gradients observed in many 
clusters. These works show that when radiative cooling is 
also included, the models with and without preheating pro- 
duce nearly indistinguishible radial entropy profiles. In all 
these calculations, regardless of the treatment of cooling, 
the large isentropic cores vanish. It is interesting to note 
that in all cases the non-preheated scenarios provide the 



most satistifactory description of the our XMM-Newton 
data, predicting entropy distributions with slopes close to 
our observed of 0.92+q q5 within O.Srsoo- In contrast, de- 
pending on the redshift at which the energy injection oc- 
curs, the preheated models predict shallower inner slopes 
than we observe. Thus the "similarity break" in the cluster 
luminosity-temperature relation — one of the original mo- 
tivations for the notion of preheating — continues to elude 
self-consistent modeling. 

The second cliici disagreement — the entropy deficit at 
the outer radii near r^oo — is more difficult to explain. 
While it is possible that substructure or superposition with 
background clusters could cause the deviation from a sin- 
gle power law, we do observe the decrement in groups 
that show no evidence of irregularity, such as SS2bl53, 
NRGbl84 and NGC 1550 (Sun et al. 2003). An under- 
standing of this entropy deficit will likely require higher 
signal-to-noise data as well as preheating models that more 
accurately reflect the temperature distribution of the in- 
tragroup medium. 

0.25 
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Fig. 8. — The pressure transition slope Tp eompared to the entropy 
transition slope rs- The solid straight line corresponds to rs = rp\ 
cirrows indicate upper or lower limits. 



Returning to the self-similar fits, we briefiy comment 
on the physical significance of the independently derived 
transition radii and Tp. Because these breaking radii are 
typically » 10 kpc, it is unlikely that they are simply re- 
lated to the effect of the potential of the central galaxy on 
the gas distribution. Similarly, most of the measured radii 
are larger than the size of a typical group-scale cooling 
flow (« 20 kpc). Thus, it is more likely that the transition 
radii are reflective of the properties of the gas distribu- 
tion throughout the group, rather than just the effect of a 
central galaxy or cooling flow. In that case, given a per- 
fectly self-similar sample, one would expect a monotonic 
functional relationship between Vs and rpi as the entropy 
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Fig. 9. — (Left) Histogram of velocities for galaxies within 12' of the center of tlie RGH 80 group. Tliere is a gap of 490 km s^^ between 
the two peaks, comparable to the group's overall velocity dispersion. (Right) POSS image of RGH80, showing galaxies belonging to the 
histogram's low velocity peali (squares) and high velocity peals (circles). The circles are concentrated around the SW X-ray emission peals. 




Fig. 10. — ii-band Gemini GMOS acquisition images of RGH80, showing the dominant X-ray peak (left) and the subgroup to the SW 
(right). Also shown are lines of constant emission-weighted temperature. In each case, the innermost contour shows a temperatmre of 0.6 
keV; each larger contour corresponds to an increment of 0.05 keV from the previous one. 
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transition radius increases or decreases, the pressure ra- 
dius should follow suit. The simplest such scenario would 
be = Tp, i.e., the entropy and the pressure profiles both 
breaking at the same point. We examine this possibility in 
Figure 8. Here we sec that half the sample docs lie close to 
the rs = rp relationship: RGH 80, NGC 5129, SS2bl53, 
and Abell 194 are within 1 cr of the line. Of these, the 
first three appear relaxed in both their X-ray and optical 
morphology (see §5 below). Of systems that deviate from 
the Ts = r-p line, Abcll 2634, SRGbll9, and HCG 97 show 
hints of irregularities that we discuss in §5 below, while 
NRGbl84 appears relaxed. Thus the relationship between 
rs and rp is a useful (if not foolproof) tool in evaluating 
the equilibrium state of a system of galaxies. 

Thus the XMM-Ncwton observations demonstrate that 
our sample of groups has a self-similar azimuthally aver- 
aged pressure distribution. The entropy profiles also ap- 
pear similar. These entropy profiles have shapes that chal- 
lenge detailed theoretical models that include preheating 
(Tozzi et al. 2000; Tozzi & Norman 2001). Furthermore, 
non-preheated models are valid only in the inner regions — 
we find S{r) oc 7-0.9i±o.04^ closely matching the the models' 
predicted S{r) oc r^'^. However, as our data approaches 
rsoo, the slope flattens, even in groups where the gas dis- 
tribution appears regular. Present preheating models do 
not account for the 7^ w 0.4 slope we find near rsoo- 

5. PROPERTIES OF INDIVIDUAL SYSTEMS 

Before we discuss the properties of the individual sys- 
tems, the heterogeneous nature of the sample as a whole 
is worth pointing out. No two systems in our sample; re- 
semble one another morphologically. There is a spherical 
and relaxed-seeming group (SS2bl53), a Hickson (1982) 
compact group merging with a sub-group (HCG 97), an 
apparent merger that the optical data reveal to be a super- 
position (RGH 80 and Abell 194), a hnear cluster (Abell 
194), and three loose groups. These systems are a random 
and plausibly representative subsample of the optically se- 
lected RASSCALS galaxy groups. 

5.1. RGH80 

This group of galaxies was originally described by 
Ramella et al. (1989); it was observed by ROSAT as well 
as ASCA, which yielded spectra with emission-weighted 
temperatures 1.5 keV and metallicitics » 0.5 solar 
(Hwang et al. 1999). The XMM observations of this group 
are described in detail by Xue, Bohringer, & Matsushita 
(2004), who find that the data are consistent with a mono- 
tonically decreasing metallicity distribution. 

Our new optical data suggests the presence of a second 
group in the field of RGH 80. The X-ray peak to the 
SW of the central peak contains another concentration of 
galaxies. These galax;ies arc also distinct from the main 
peak. Figure 9 shows a histogram of the rest-frame veloci- 
ties of the galaxies within the central 12'. The 21 galaxies 
in this region exhibit a 490 ± 10 km s~^ void in velocity 
space. The centers of the peaks themselves are separated 
by 1080 ± 103 km s"-'^. If we split the two systems at the 
midpoint of this void, we see that 7 of the 8 galaxies in 
the higher-velocity peak reside within 30" (20 kpc) of the 
SW X-ray peak. This lends support to the idea that the 
SW peak, separated from the main peak by 156 ± 5 kpc, is 



a dynamically distinct group of galaxies. The two groups 
are shown in Figure 10. 

Several lines of reasoning suggests that the SW group 
has not interacted with the main group. The void in the 
velocity histogram (490 km s^^) is comparable to the ve- 
locity dispersion of the group overall (~ 600 km s~^). The 
size of this void argues against a scenario in which the two 
groups have already interacted. If one group were bound 
to and had already passed through the other, at least a 
few galaxies within our detection limit ought to have been 
displaced to occupy the velocity gap. 

We can make this argument more rigorous by examining 
a two body model first used by Beers. Geller, & Huchra 
(1982) to study rich clusters of galaxies. In this simplified 
scenario, the two groups are approximately point masses 
in an orbit with zero angular momentum. Given a line-of- 
sight velocity diff'erence and a projected distance between 
the two bodies, the model allows one to solve for M{i), 
the total mass required to bind the system as a function 
of the inclination angle i between the observer's line of 
sight and the orbital axis. We plot this function in Fig- 
ure 11. We also show the mass inferred from our tem- 
perature measiiremcnts via the X-ray mass-temperature 
relation (Finoguenov, Reiprich, & Bohringer 2001). In the 
case of RGH 80, the main NE clump has an emission- 
weighted temperature in the range 0.7-0.9 keV; for the 
SW clump, the range is 0.5-0.7 keV, corresponding to 
an estimated mass mass of 1.8-3.2 xlO'^'''Af0 within r^oo- 
(68% confidence, taking into account errors in the mass- 
temperature relation as well as our own measurement er- 
rors) . To estimate the remaining material within the virial 
radius r2oo ; we assume that the matter distribution of both 
clumps can be described by a Navarro et al. (1997) pro- 
file with transition radius r^oo, we arrived at this value by 
examining the mass profiles derived in the infall regions 
of similar systems by Rines et al. (2003), and converting 
their best-fit NFW profiles to our r5oo reference frame. 
The virial mass of the two clumps is then 1.3 times the 
mass at r^oo, or 2.3-4.2 xlO^^M©. The minimum mass 
required to bind the group is about 3 times larger than 
the X-ray mass inferred from the mass-temperature rela- 
tion. In the context of the simple two-body model, then, 
it is unlikely that the groups are bound. 

The X-ray morphology leads to a similar conclusion. 
The maps of this region show two distinct entropy min- 
ima, with no shocks indicating a merging history. Both 
the central clump and the southern clump are embedded in 
separate low-entropy zones. It is more likely that the SW 
pressure enhancement corresponds to a separate dark mat- 
ter concentration associated with the potential of the SW 
group, and less likely that it corresponds to a shock. The 
spectral fits show that there are two components consis- 
tent with being projected on top of each other. In the NE 
source, we observe a rise of the temperature from 0.9 keV 
to 1.2 keV, and then in the outer regions back down again 
to 0.7 keV. The metallicity is enhanced in the SW peak 
compared to the rest of the image — it is distinct from its 
surroundings. Thus the X-ray observations also support a 
scenario in which the SW group is dynamically indepen- 
dent of the rest of the system. 
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Fig. 11. — Simple two-body analysis of tlic NE and SW groups 
in RGH 80. The two solid lines indicate upper and lower limits on 
the minimum total mass required to bind the two groups to each 
other. For inclinations i < 89, only incoming bound solutions are 
possible; for i > 89°, both incoming and outgoing bound solutions 
are possible, but the latter occupy only a miniscule region in the 
solution space. The lower shaded region indicates the total mass of 
the system derived from the X-ray mass-tcmperature relation. The 
X-ray mass is too small to bind the groups by about a fa<;tor of 5. 

5.2. Ahell 194 

Abell 194 was originally classified as a "linear" cluster 
of galaxies by Rood & Sastry (1971) because of its NE-SW 
alignment. Our XMM-Newton observation shows that the 
alignment also exists in the X-ray. The first impression is 
that there are three peaks: one to NE, one in the center, 
and one to the SW. The X-ray spectrum of the central 
15"of the NE peak is consistent with that of an active 
galactic nucleus (AGN). 

It is however the southwesternmost X-ray peak that 
would seem to be of particular interest, because Knezek 
& Bregman (1998) report that the galaxy at its center 
has a redshift consistent with that of Abell 194. Yet this 
galaxy is at least four magnitudes fainter than the galax- 
ies occupying the two other, bright X-ray emission peaks. 
This would imply an unreasonably large mass-to-light ra- 
tio of w 1000 for the SW X-ray peak. Upon reanalyzing 
the Knczck & Bregman (1998) spectrum (Figure 12), wc 
find that its recession velocity is 47900±140 km s^i. Also, 
the X-ray spectrum corresponding to the SW region is bet- 
ter fit by a thermal plasma at redshift 0.15 than it is by a 
plasma at redshift 0.017. The galaxy atop the SW X-ray 
peak is therefore the brightest member of a more distant 
cluster, and the peak is not associated with Abell 194, 
but corresponds to the X-ray atmosphere of the cluster at 
z = 0.15. This result, together with the result for RGH 
80, stresses the importance of projection effects when con- 
sidering the properties of clusters and groups of galaxies 
(Gen 1997). 

Thus, there are only two distinct halos in Abell 194, 

both exhibiting an outward rise in temperature from 0.7 
keV in the center to 1.6 keV in the outskirts. The velocity 
distributions of the galaxies within each X-ray peak are in- 
distinguishable; as a result it is likely that the two clumps 
are bound. 



The entropy profile of Abell 194 seems to resist an az- 
imuthally averaged description (§ 4). The cluster registers 
a larger than the other seven systems in our sample, 
regardless of the entropy model. Given that it is by far 
the least axisymmctric cluster in our sample, the disagree- 
ment is hardly a surprise. However, it is surprising that 
the pressure profile is fully consistent with the self-similar 
model fit to the rest of the sample. 

5.3. HCG 97 

HGG 97 is a group of galaxies originally discussed in 
the Hickson (1982) catalog of compact systems. Ebeling 
et al. (1994) first reported X-ray emission in the ROSAT 
All-Sky Survey for this group. 

Our analysis of the XMM-Newton observation reveals a 
complex structure in this interacting system, with a large 
gas tail stretching toward the SE. This tail is a source of 
low-entropy gas and may account for the deficit seen in 
the entropy profile at around 0.5r5oo- The Gemini obser- 
vation reveals a single member galaxy inhabiting the SE 
extension plume; a close examination of the GMOS acqui- 
sition image (Figure 13) reveals no background clusters 
that could be responsible for the X-ray emission in that 
region. 

If the SE plume is associated with the cluster, it may 
be be a relic of the infall of one of the brighter interact- 
ing galaxies at the group center — or it may be yet another 

chance superposition from a cluster not visible in the Gem- 
ini image. The temperature of the plume is consistent with 
that of its surroundings: i.e., the plume is not apparent in 
the temperature map (Figure 4). However, the plume is 
clearly present in the entropy, pressure, and density dis- 
tributions. Neither the X-ray and the optical data are at 
this point of sufficient quality to determine whether the 
plume has been stripped from one of the central galaxies, 
or whether it is a background structure unrelated to the 
main group. 

Substructure notwithstanding, the azimuthally aver- 
aged pressure and entropy distributions in HGG 97 are 
consistent with the self-similar power laws fitted to the 
sample as a whole. 

5.4. SS2bl53 

SS2bl53 is a remarkably round system, a stark con- 
trast to the complexity of the other groups. Because 
the only redshifts available in its region are from the 
nearly decade-old Las Gampahas Redshift Survey (Sheet- 
man et al. 1996), very few galaxies are known to be mem- 
bers. 

The azimuthally averaged pressure and entropy distri- 
butions in this system are consistent with the self-similar 
broken power law fit to the sample as a whole, except for 
the innermost data point. The pressure at w 10 kpc dis- 
agrees with the self-similar law by a factor of 2. SS2bl53 
is the most nearby group in the sample, and it is possible 
that the pressure enhancement represents the contribu- 
tion of the gas bound to the central galaxy, NGG 3411, 
to the total pressure. Similar enhancements — associated 
with the central galaxies of groups within z « 0.01 — have 
been observed in other systems (Helsdon et al. 2001). 
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al94gal.fits]l comb vel= 47985.94 (52.96) km/sec R= 9.57 
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Fig. 12. — Continuum-subtracted spectrum of the galaxy associated with the SW emission peak in the XMM-Newton observation of Abell 
194. The redshift of the galaxy indicates that it is a distant background cluster, and not related to Abell 194. 




Fig. 13. — (left) A 10 minute, R band Gemini GMOS acquisition image of HCG 97. X-ray surface brightness contours from our wavelet 
decomposition technique are superimposed. The plume to the SE {right) is coincident with a single member galaxy. There axe no noticeable 
background clusters associated with the peak within the SE extension. 
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5.5. NGC 5129 Group 

The "plume" to the NE seems to be indistinguishable 
from its environment to within the errors. The entire cen- 
tral region seems roughly isothermal, with a temperature 
0.7 keV. The entropy profile of the NGC 5129 Group does 
not reveal any deviations from the self-similar models, and 
its pressure profile also appears consistent. 



srgbll9.034.flts]l comb vel= 5570.10 (29.95) km/sec R= 19.80 



5.6. Ahell 

The XMM-Newton observation of Abell 2634 was 
plagued by long periods of high X-ray background, which 
did not allow us to extract as much spectroscopic informa- 
tion as would have otherwise been possible for this bright 
cluster. The radio source 3C 465 coincides with the cen- 
tral galaxy (Owen & Rudnick 1976; Sakelliou & Merrifield 
1999); it exhibits a powerful jet that interacts the X-ray 
emitting gas as far out as 100 kpc from the cluster center. 
For this reason, our entropy and pressure profiles miss the 
central section of the cluster entirely. 



hut 



_L 
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srgbll9.035.fitsjl comb vel= 5984.78 (36.60) km/sec K= 14.53 
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Fig. 14. — Continuum-subtracted spectra of the western (top) 
and eastern (bottom) pair of central galaxies in SRGbllQ. . 



5.7. NRGbl84 

The group has a north/south elongation in surface 
brightness. This elongation is less pronounced in the pres- 
sure map. There is a ~ 20% pressure difference between 
the regions directly to the SE and NW of the group center. 
The temperature rises fairly quickly from 0.6 keV in the 
very central region to about 1.5 keV. The N-S extension 
that is visible in the surface brightness map is not present 
in temperature map. On both the entropy and pressure 
plots, the properties of the gas in NRGbl84 appear typical. 



5.8. SRGhllQ 

In general, the temperature increases from 0.7 keV to- 
wards the center to 1.4 keV towards the exterior. There 

are too few counts to constrain the iron abundance ro- 
bustly. - There is an interesting unresolved source to the 
east of the very center of the group. It is a member galaxy 
that is ~ 2 mag fainter than the brightest group galaxy. 
The two galaxies have a large velocity difference (450 km 
s~^). Despite being fainter, the galaxy to the E has a 
higher temperature than the central source, 1.6 ± 0.2 keV 
compared to 0.7 ± 0.2 keV. The optical spectra of both 
the central galaxies show no evidence of AGN emission 
(Figure 14). The higher X-ray temperature for the source 
to the east is therefore a puzzle. Our analysis rules out a 
significant contribution coming from low-mass X-ray bina- 
ries (LMXBs): the spectrum is inconsistent with the > 5 
keV tempcratmes expected from these objects (Sarazin, 
Irwin, & Bregman 2001). At the same time, however, a 
power law spectrum with index 2.9 ± 0.4 provides a better 
fit than a thermal plasma. The power law spectrum has a 
total luminosity of 9 x 10^° ergs/s in the 0.5-2 keV band. 
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The source of this spectrum could be a veiled AGN such 
as the one observed in NGC 1291 by Irwin, Sarazin, & 
Bregman (2002). 

The pressure and entropy distribution in SRGbllQ 
presents a puzzle similar to that of AbcU 194. By itself, 
the azimuthally averaged entropy profile of SRGbll9 is 
not a good fit to the simple empirical models in §4. Of the 
eight groups in our sample, SRGbll9 makes the second 
largest contribution to the reduced (after Abell 194), 
regardless of the shape of the model. However, similarly 
to Abell 194, SRGbll9 is quite consistent with the self- 
similar pressure models. Thus we find that the two groups 
with the more discrepant entropy distributions are never- 
theless similar to the rest of the sample when we consider 
the pressure profiles. 

6. CONCLUSION 

We use new and previously unpublished archival XMM- 
Newton data to examine the properties of the X-ray emit- 
ting medium in eight groups of galaxies from the RASS- 
CALS survey. In addition to the X-ray data, we collect 
redshifts for the group members from the literature. We 
also report redshifts from a Gemini North Observatory 
survey of the groups RGB 80, BCG 97, and NRGbl84. 

The sample consists of groups with X-ray luminosities 
< 10^^ erg s~^, and was chosen based on availability in the 
XMM-Newton data archive. The groups in the sample ex- 
hibit a wide variety of X-ray and optical morphologies, 
from seemingly perfectly spherical systems (SS2bl53) to 
linear clusters (Abell 194). Many of the systems exhibit 
certain or possible superposition, in both the X-ray and 
optical, with unrelated background galaxies or clusters of 
galaxies. We use the Gemini North optical data to ar- 
gue that RGB 80 is a superposition of two dynamically 
unrelated groups. 



Despite the morphological incongruity of the groups, 
they exhibit regularity in their azimuthally averaged pres- 
sure and entropy profiles. Self-similar broken power laws 
can describe the spatial distribution of both these quanti- 
ties. For the entropy, which increases with distance from 
the center, we find an inner slope = 0.921o;q5 and outer 

slope 7, = 0A2trol (68% confidence). These results ar- 
gue against the existence of an iscntropic "floor" at the 
center of X-ray emitting groups, and against theoretical 
models that describe the entropy distribution as a single 
power law for r > O.Srsoo- 

The decreasing pressure profiles are also well-fitted by 
a self-similar broken power law. The inner and the outer 
slopes are ap = — 0.78lo!o3 and jp = — I.TIq^s 
tire sample. The regularity and similarity of the pressure 
profiles suggests that the groups reside in dark matter ha- 
los that are similar in structure. 

We have made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Admin- 
istration. We thank Patricia Knezeck for providing us 
with the spectrum of the galaxy in the SW X-ray peak in 
Abell 194, and the anonymous referee for comments that 
improved the paper. J. P. H. thanks the Alexander von 
Bumboldt Foimdation for support in Garching. M. J. G. 
was supported by the Smithsonian Institution. A. F. ac- 
knowledges support from BMBF/DLR under grant 50 OR 
0207, from MPG, and from NASA grant NNG04GB226. 
A. M. was supported by NASA through a Chandra Post- 
doctoral Fellowship Award issued by the Chandra X-ray 
Observatory Center, which is operated by the Smithsonian 
Astrophysical Observatory for and on behalf of NASA un- 
der contract NAS 8-39073. 
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Table 1 
Group Sample 



Primary ID 02000 ^2000 AT cz crios Other ID Central Galaxy Sources 



Abell 194 


01:25:48.1 


-01:23:06 


106 


5316t| 


550111 


SRGbl03 


a 


R(106) 


SRGbll9 


01:56:20.8 


+05:37:40 


61 


5525_^g 




WBL 61 


NGC 741 


M(21), Z{30) 


SS2bl53 


10:50:26.3 


-12:50:43 


5 


4747+1^ 


16lt« 


USGC 8152 


NGC 3411 


L{14) 


NRGbl84 


12:08:05.5 


-1-25:14:14 


33 


6700l« 


390t^? 




UGC 07115 


G(8), M(52) 


RGH80 


13:20:14.8 


-1-33:08:35 


58 


110471^^ 


602+^f 


NRGs241 


NGC 5098 


G(14), M(44) 


NGC 5129 


13:24:09.6 


-1-13:58:51 


33 


69521^6 


283+2? 


NRGb244 


NGC 5129 


M{19), Z(15) 


AbcU 2634 


23:38:29.2 


-1-27:01:54 


69 


939ll72 


5681^9 


3C 465 / SRGs040 


NGC 7720 


P(69) 


HCG 97 


23:47:23.1 


-02:18:07 


37 


6758lg 


383t|° 


SS2b312 


IC 5357 


G(ll), M(26) 



Note. — Af is the number of known members within 2 Mpc of the X-ray center of the group; cz is the mean redshift of the group 
multiplied by the speed of light in km s~^; a\o^ is the velocity dispersion of the group in km s~^. Also shown are other names of 
each group, as well as the common name of the central galaxy. The redshifts are taken from the following sources: Mahdavi & 
Geller (2004, M); Shectman et al. (1996, L); Zabludoff & Mulchaey (1998, Z); Ramella et al. (2002, U); Rines et al. (2003, R); 
Pinkney et al. (1993, P); and Gemini (G) redshifts measured in this work. The number after each letter indicates the number of 
unique galaxies that were measured in that work. 

■^The "chain cluster" Abell 194 does not have a dominant galaxy. 
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Table 2 

Xew" Optical Data 





(52000 


cz 




TvDe 


12:07:44.31 


+25:14:48.5 


41241 


8 


E 


12:07:45.67 


+25:14:36.9 


26338 


5 


E 


12:07:46.99 


+25:15:51.5 


6468 


16 


A 


12:07:47.71 


+25:11:13.2 


62187 


102 


A 


12:07:49.53 


+25:12:35.9 


29447 


9 


E 


12:07:54.08 


+25:16:58.9 


155457 


22 


E 


12:07:55.33 


+25:14:00.2 


89416 


30 


A 


12:07:57.64 


+25:16:43.3 


126564 


12 


E 


12:07:58.16 


+25:10:17.8 


112008 


38 


A 


12:07:59.33 


+25:11:45.9 


5896 


59 


A 


12:08:00.12 


+25:15:14.9 


85701 


25 


E 


12:08:01.94 


+25:15:43.3 


100206 


40 


A 


12:08:02.23 


+25:11:18.2 


7096 


36 


A 


12:08:04.01 


+25:16:30.1 


134066 


36 


A 


12:08:04.68 


+25:19:10.0 


142108 


24 


E 


12:08:06.12 


+25:16:02.9 


5813 


23 


A 


12:08:06.23 


+25:13:31.6 


160608 


44 


A 


12:08:07.16 


+25:15:36.9 


6502 


19 


A 


12:08:07.44 


+25:12:51.5 


6964 


38 


A 


12:08:08.16 


+25:14:44.3 


89252 


22 


A 


12:08:08.86 


+25:13:05.2 


25449 


19 


A 


12:08:10.42 


+25:16:03.3 


197516 


30 


E 


12:08:11.97 


+25:11:21.9 


24130 


12 


E 


12:08:12.37 


+25:18:04.3 


189232 


42 


E 


12:08:13.18 


+25:15:40.3 


92668 


46 


A 


12:08:14.87 


+25:09:38.8 


153102 


27 


E 


12:08:15.42 


+25:13:23.3 


155169 


31 


E 


12:08:15.86 


+25:16:18.6 


79235 


14 


E 


12:08:17.88 


+25:14:32.2 


78955 


16 


E 


12:08:18.85 


+25:13:46.3 


78683 


9 


E 


12:08:19.19 


+25:15:47.2 


189886 


34 


E 


12:08:19.37 


+25:13:57.9 


6912 


60 


A 


12:08:21.44 


+25:09:40.9 


110299 


14 


E 


12:08:22.19 


+25:17:58.6 


25738 


13 


E 


12:08:22.26 


+25:18:57.8 


6827 


26 


A 


12:08:22.88 


+25:09:48.9 


110235 


17 


E 


12:08:25. .33 


+25:16:20.9 


78222 


15 


E 


13:20:00.91 


+33:09:11.9 


129328 


61 


A 


13:20:07.44 


+33:07:00.8 


10991 


60 


A 


13:20:08.12 


+33:10:38.0 


10999 


29 


A 


13:20:08.32 


+33:06:13.0 


12317 


55 


A 


13:20:14.00 


+33:07:34.0 


10483 


51 


A 


13:20:14.70 


+33:12:17.6 


10594 


27 


A 


13:20:15.01 


+33:05:38.1 


11922 


6 


E 


13:20:15.05 


+33:11:57.3 


67931 


28 


E 


13:20:15.99 


+33:04:38.3 


29945 


11 


E 


13:20:17.45 


+33:08:54.4 


80238 


81 


A 


13:20:17.76 


+33:12:12.1 


35837 


10 


E 
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Table 2 — Continued 



0^2000 


02000 


cz 




Type' 


13:20:18.79 


+33:10:17.9 


105272 


67 


A 


13:20:18.80 


+33:09:55.0 


104821 


12 


E 


13:20:19.03 


+33:08:03.9 


10725 


33 


A 


13:20:20.01 


+33:13:30.1 


11005 


14 


E 


13:20:21.27 


+33:04:41.0 


10283 


6 


E 


13:20:23.04 


+33:09:17.5 


11001 


65 


A 


13:20:23.59 


+33:06:40.3 


98242 


34 


A 


13:20:25.40 


+33:11:24.2 


5460 


20 


E 


13:20:27.84 


+33:08:27.4 


10735 


27 


A 


13:20:31.15 


+33:06:08.8 


68601 


29 


A 


13:20:31.69 


+33:07:14.3 


11008 


62 


A 


13:20:33.30 


+33:05:59.0 


11714 


5 


E 


13:20:33.40 


+33:12:00.0 


141341 


29 


E 


13:20:35.19 


+33:05:40.7 


133709 


41 


E 


13:20:36.30 


+33:09:45.0 


20144 


17 


E 


23:47:05.53 


-02:14:30.6 


42715 


4 


E 


23:47:06.20 


-02:14:10.4 


35609 


21 


E 


23:47:09.34 


-02:15:37.6 


7568 


48 


A 


23:47:09.37 


-02:22:40.8 


89460 


22 


E 


23:47:09.92 


-02:18:54.8 


140450 


18 


E 


23:47:10.10 


-02:15:11.8 


7591 


6 


E 


23:47:10.65 


-02:19:38.7 


116770 


23 


E 


23:47:11.65 


-02:22:34.3 


152617 


46 


A 


23:47:14.70 


-02:21:53.4 


45685 


14 


E 


23:47:15.40 


-02:21:35.5 


115631 


54 


A 


23:47:15.60 


-02:20:50.4 


115784 


25 


A 


23:47:16.20 


-02:19:09.0 


96754 


7 


E 


23:47:18.90 


-02:24:12.8 


77719 


10 


E 


23:47:19.40 


-02:21:19.0 


89411 


26 


E 


23:47:19.90 


-02:16:50.4 


6665 


24 


A 


23:47:19.94 


-02:16:34.0 


63116 


14 


E 


23:47:20.10 


-02:18:37.6 


88188 


25 


A 


23:47:20.30 


-02:22:26.8 


7247 


4 


E 


23:47:20.94 


-02:17:32.5 


124893 


24 


E 


23:47:21.00 


-02:19:28.5 


140106 


32 


E 


23:47:21.70 


-02:17:32.2 


29293 


11 


E 


23:47:22.00 


-02:18:53.4 


29389 


53 


A 


23:47:22.00 


-02:21:12.6 


6033 


48 


A 


23:47:22.62 


-02:18:14.6 


68783 


39 


A 


23:47:22.80 


-02:15:25.0 


88275 


30 


A 


23:47:23.30 


-02:21:47.0 


5991 


7 


E 


23:47:23.88 


-02:16:02.6 


111777 


28 


E 


23:47:24.72 


-02:16:35.6 


152587 


14 


E 


23:47:24.96 


-02:16:24.1 


124710 


38 


A 


23:47:25.00 


-02:19:35.5 


5837 


37 


A 


23:47:26.32 


-02:16:44.4 


55222 


8 


E 


23:47:26.60 


-02:21:11.7 


6778 


36 


A 


23:47:26.60 


-02:21:11.8 


6519 


25 


E 



Table 2 — Continued 



0:2000 ^2000 cz ecz Type* 



23:47:26.60 


-02: 


21:11.8 


6698 


6 


E 


23:47:27.36 


-02: 


17:22.9 


169001 


32 


E 


23:47:28.36 


-02: 


24:00.3 


169622 


27 


A 


23:47:28.81 


-02: 


20:12.7 


-117 


45 


A 


23:47:29.50 


-02: 


18:04.8 


68533 


16 


E 


23:47:29.50 


-02: 


18:05.0 


68523 


12 


E 


23:47:29.60 


-02: 


22:56.3 


169789 


5 


E 


23:47:30.00 


-02: 


23:30.8 


41439 


13 


E 


23:47:30.23 


-02: 


10:51.1 


125224 


26 


A 


23:47:30.50 


-02: 


22:01.6 


168233 


30 


E 


23:47:30.80 


-02; 


;19;08.3 


45611 


30 


A 


23:47:31.20 


-02; 


; 20:28. 2 


42185 


8 


E 


23:47:31.20 


-02: 


; 20:28.2 


42206 


45 


A 


23:47:32.10 


-02: 


;18:18.3 


72077 


16 


E 


23:47:32.40 


-02: 


;22:29.3 


6879 


66 


A 


23:47:32.47 


-02: 


;16:24.8 


115977 


39 


A 


23:47:34.65 


-02; 


;18:48.5 


185245 


28 


E 


23:47:35.20 


-02; 


;21:39.1 


84656 


17 


E 


23:47:35.74 


-02; 


;15:30.1 


135273 


37 


E 


23:47:36.00 


-02; 


;24:10.7 


129738 


26 


A 


23:47:36.83 


-02; 


;17:13.0 


69053 


9 


E 


23:47:38.20 


-02; 


;20:38.0 


183269 


24 


E 


23:47:38.47 


-02; 


:16:28.4 


106316 


28 


E 


23:47:39.27 


-02; 


:14:38.2 


169012 


10 


E 


23:47:40.04 


-02; 


18:57.9 


139059 


11 


E 


23:47:40.90 


-02; 


21:01.5 


68338 


12 


E 


23:47:42.40 


-02; 


:23:47.5 


66326 


19 


E 


23:47:42.64 


-02; 


14:33.1 


115187 


21 


A 


23:17:42.70 


-02; 


:21:ol.2 


129539 


12 


A 


2:-!: 17:1:-!. 80 


-02; 


17; 10. (i 


l:-!S877 


21 


E 



Note. — cz denotes the redshift of the object multi- 
phed by the speed of Ught; tcz is the la uneertainty in 
cz; "type" indicates whether the dominant features in 
the spectrum are absorption ("A") or emission ("E") 
lines. 



Table 3 

XMM Observational Information 



Name 


Obs. ID 


net exp filter 


Orbit 


frame time 






ksec 






ms 


AbcU 194 


0136340101 


12.9 


Thin 


557 


73 


SRGbllQ 


0153030701 


3.9 


Thin 


745 


73 


SS2bl53 


0146510301 


18.5 


Thin 


555 


73 


NRGbl84 


0151400201 


8.6 


Thin 


634 


73 


RGH 80 


0105860101 


26.2 


Thin 


563 


199 


NGC 5129 


0108860701 


12.4 


Medium 


472 


73 


Abell 2634 


0002960101 


4.0 


Medium 


464 


73 


HOG 97 


0152860101 


20.9 


Thin 


559 


73 



23 



Table 4 



IGM Properties of the Groups between 0.1 — O.Srsoo- 



Group 


''500 


{kT) 


(Z) 


(S) 


{P) 




kpc 


keV 


Zq 


keV cm^ 


10 ergs cm ^ 


Abell 194 


521 


1.36 ±0.04 


0.19 ±0.03 


260. ± 15.1 


0.89 ±0.06 


SRGblig 


520 


1.36 ±0.07 


0.55 ±0.27 


402. ± 98.7 


0.46 ±0.13 


SS2bl53 


407 


0.65 ±0.02 


0.15 ±0.01 


109. ± 6.8 


0.47 ±0.03 


NRGbl84 


488 


1.26 ±0.07 


0.07 ±0.02 


192. ± 13.0 


1.11 ±0.09 


RGH 80 


455 


1.04 ±0.02 


0.16 ±0.02 


107. ±6.1 


1.67±0.10 


NGC 5129 


434 


0.95 ±0.05 


0.39 ±0.12 


178. ±26.2 


0.60 ±0.13 


Abell 2634 


639 


2.07 ±0.11 


0.20 ±0.05 


259. ±16.3 


2.42 ±0.16 


HOG 97 


420 


0.89 ±0.02 


0.17 ±0.03 


129. ±8.6 


0.79 ±0.06 



Table 5 

Self-similar Pressure and Entropy Profile Fits 

Entropy Fits Pressure Fit 
Entire Sample 

Goodness-of-fit 0.0095 (x^/i^ = 93/62) 0.15 {x^/f = 73/62) 

as 0.87-0.96 -ap 0.75-0.82 

7s 0.38-0.47 -7p 1.4-1.8 
Abell 194 

Ts >0.17 rp 0.16-0.20 

So 120-140 Po 17-20 
SRGbll9 

rs >0.14 rp 0.042-0.056 

So 150-190 Po 17-26 
SS2bl53 

rs 0.054 0.069 rp 0.014 0.047 

So 100-120 Po 19-23 
NRGbl84 

Vs 0.062-0.083 rp > 0.19 

So 170 200 Po 12-15 
NGC 5129 

Vs 0.062-0.10 rp 0.042-0.070 

So 160-200 Po 11-14 

RGH 80 

Vs 0.048-0.068 rp 0.042-0.096 

So 110-130 Po 26-31 
Abell 2634 

rs < 0.062 rp 0.12-0.14 

So 240-600 Po 75-96 
HCG 97 

rs 0.021-0.030 rp 0.083-0.098 

So 180-250 Po 14-16 



Note. — the 68% confidence intervals take into account correlations 

with all the other parameters. The transition radii and rp are in units 
of Mpc; 5*0 is in units of keV cm^, and Pq is in units of 10"^^ erg cm^'^. 



